Low-frequency current fluctuations in doped ladders 
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Charge current static and dynamical correlations are computed by exact diagonalisation methods 
on a 2-leg t-t'-J ladder which exhibits a sharp transition between a Luther-Emery (LE) phase of hole 
pairs and a phase with deconfined holes. In the LE phase, we find short-range low-energy incommen- 
surate current fluctuations which are intrinsically connected to the internal charge dynamics within 
one hole pair. On the contrary, when holes unbind, the maximum of the current susceptibility moves 
abruptly to the commensurate wavevector ir and strongly increases for decreasing doping suggesting 
an instability towards a staggered flux state at sufficiently small doping. 

PACS: 75.I0.-b, 75.50.Ee, 71.27. +a, 75.40.Mg 
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Among mean-field solutions for two-dimensional (2D) 
strongly correlated electrons, the staggered flux phase ||] 
appeared as an appealing candidate. Due to the local 
SU(2) symmetry, the staggered flux phase is equivalent, 
at half-filling, to the projected d-wave pairing BCS varia- 
tional wavefunction Q . Both states could equivalently be 
viewed as a Resonating Valence Bond state || . Since hole 
doping lifts the degeneracy between the staggered flux 
state and the d-wave pairing state, it is of importance to 
investigate possible competition between such states. At 
finite doping, as the newly proposed d-density wave state 
j|, the staggered flux state would break translation as 
well as time reversal symmetry. In otherwords, it would 
exhibit "orbital antiferromagnetism" i.e. a checkerboard 
pattern of plaquette currents circulating clockwise and 
anti-clockwise on the two different sublattices. On the 
other hand, it has been shown ||] recently that the cur- 
rent vorticity correlations in a variational d-wave pairing 
state follows closely the hole correlations within one pair. 

Numerical investigations carried out on small (2D) 
clusters have shown that staggered or commensurate 
(with a larger period modulation) flux states are very 
good variational states pj. In addition, spin-chirality 
correlations have been tested numerically in the 2D t- 
J model @ showing that staggered (i.e. at q = (7r,7r)) 
chiral fluctuations are strongly enhanced with increas- 
ing hole kinetic energy. Recent calculations of the static 
current-current correlations on a 32-site lattice with two 
holes H have confirmed the existence of staggered vor- 
ticity in this model. 

In this paper, search for current fluctuations have been 
carried out in an extended t-J model on a two-leg lad- 
der where finite size effects can be handled in a controled 
fashion. Such a system is known to exhibit a robust spin 
gap at and close to half-filling as well as hole pair bind- 
ing Dominant power-law d x 2_ a 2-like pairing and 
Akp charge density wave (CDW) correlations at small 
doping are characteristic of a Luther-Emery (LE) liquid 
regime JTl[ . An additional negative diagonal hopping |l2] ] 
t' between next-nearest neighbor (NNN) sites was shown 
to strongly reduce hole pair binding leading e.g. to the 



emergence of low-energy quasiparticle-like triplet excita- 
tions . On the contrary, a positive t' increases the pair 
binding energy even further. Here, we use this parameter 
to monitor the hole pair binding energy in order to inves- 
tigate possible connections between staggered-vorticity 
correlations and hole correlations. We find that incom- 
mensurate current fluctuations occur on the scale of the 
pair size. For sufficiently negative t' hole pairs break ap- 
part |l| and low energy staggered vorticity fluctuations 
are found as well as static commensurate correlations in 
the ground state (GS). We argue that these finding might 
signal an instability towards a staggered flux state. 
The two-leg ladder Hamiltonian reads, 
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where i is a site index along the chain direction, I and 
m (=1,2) label the 2 legs and n^i is the local charge 
density operator. Holes are assumed to hop along the 
chain direction or between legs with amplitudes t\\ +1 = 
t\ es and t\f = tfl = t rung respectively. In addition, as 
stated above, it is of particular interest to include a NNN 
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t' [pL In the following, for 



sake of simplicity, we shall only discuss the case of a 
spatially isotropic ladder i.e. with t\ cg — i rung = t (set 
to 1) and Ji e g = J rung = J. In addition, J is set to 
0.5. However, we believe our results are generic. The 
static and dynamical quantities reported in this paper 
have been obtained by exact diagonalisations (ED) of 
finite periodic 2xL clusters (L=8 to L=12) doped with 
up to 6 holes. 

The NNN hopping t' is the key parameter of our study 
since it monitors the attraction between holes in the sys- 
tem [ pL3| , To demonstrate this fact, as a preliminary 
study, we have computed the (static) charge correlations 
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on 2 hole-doped ladders for many different parameters t'. 
Selections of these data reported in Fig. [l] show clearly 
2 different behaviors; above a critical value t' c (we find 
t' c ~ —0.3 at J = 0.5) the two holes are bound and the 
hole-hole correlation decreases very fast with distance. 
Note that, in this case, finite size effects remain quite 
small since hole correlations at short distance are rather 
independent on system size. The bound pair "size" in- 
creases with decreasing t 1 and eventually at t' = t' c (neg- 
ative) the GS changes abruptly Jl4| into a state where 
holes are unbound. Indeed, as shown for t' = —0.4 in 
Fig. [j], the hole-hole correlation function increases with 
distance which means that holes tend to avoid each other. 
It is interesting to notice that boundary effects become 
quite important in this case probably in connection with 
the enhancement of the antiferromagnetic (AF) correla- 
tions. In order to avoid magnetic frustration, the 2 holes 
tend to be located, on average, on the same (opposite) 
leg for even (odd) ladder lengths 15 1. In the following, 
we shall show that the charge current fluctuations also 
change rapidly at the transition between bound pairs and 
individual holes. 



2 holes J=o.5 

2x12 (filled) & 2x11 (open) 




FIG. 1. Hole-hole (equal-time) correlation function on a 
two hole-doped 2x12 ladder vs separation between holes in 
the same leg (a) and in different legs (b). Different filled 
symbols correspond to different t' parameters according to 
the legend. Data obtained on a 2x11 ladder are also shown 
as open symbols (unless undistinguishable from the previous 
sets). Data are normalized to the hole density i.e. such that 
the local (r = 0) correlation equals 1. 

The current operators for the rung and leg bonds are 
defined as follows. For the leg, 



-3 + 1,1,^3,1,0 
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where I and n denote the site numbers for the leg and 
rung directions, respectively. Note that the leg current is 



defined in such a way that current densities on opposite 
bonds on the two legs have opposite signs as it should for 
loop currents defined on the plaquettes. For the rung, 
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In the staggered flux phase |Q both rung and leg cur- 
rents would acquire a finite expectation value at momen- 
tum q = 7r corresponding to a staggered arrangement of 
the loop currents on the plaquettes. Here we consider a 
more general situation of possible incommensurate modu- 
lations with arbitrary q flltf of the current pattern along 
the leg direction. Dynamical fluctuations of such loop 
currents can be investigated by computing the associated 
complex dynamical susceptibility, 
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where \n) (E n ) are exact eigenstates (eigenenergies) of 
the hamiltonian. It is instructive to define the current- 
current dynamical correlation function Sjj as the imagi- 
nary part of the susceptibility devided by frequency, 
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E n — Eq 



-S (u - E n + Eo) , (5) 



which leads, after integration over frequency, to the static 
(to = 0) current susceptibility (Kramers-Kronig), 
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The current-current dynamical correlation Sjj and the 
related static susceptibility have been computed on fi- 
nite clusters by a standard continued-fraction method. 
Note that the susceptibility x(q> 0) can alternatively be 
calculated as a response to a fictitious modulated mag- 
netic flux. We have checked that this second method 
gives identical results. We also have found that, quite 
generally, use of bond and leg currents gives very simi- 
lar results (especially at q = tt) so that in the following, 
for sake of clarity, only results using the rung current 
operator will be presented. 

We first focus on a ladder doped with only two holes. 
A selection of our results obtained on a 2x12 ladder are 
presented on Fig. g. On should distinguish two interest- 
ing features on these data; (a) some high energy struc- 
ture in the range w ~ it which is merely dispersing and 
whose weight increases with increasing t'\ (b) a low en- 
ergy structure at a fraction of J. Above t' c ~ —0.3, the 
maximum weight of this structure is spread on incom- 
mensurate 16| ] wavevectors in the vicinity of n. On the 
contrary, when the two holes become unbound (t' < t' c ), 
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FIG. 2. Dynamical correlation function Sjj(q,u)) calcu- 
lated on a 2 hole-doped periodic ladder for J = 0.5. Each 
panel corresponds to a given value of t as indicated. For 
clarity, curves obtained for all available momenta q are shifted 
upwards w.r.t. each other by an equal amount. From bottom 
to top, q — 2nn/L, n from to L/2. Note the reduction 
factor of 0.4 applied to the data at t 1 = —0.3 and q — n. 

the largest low energy (typically cj typ ~ 0.2 J = O.lt) 
current fluctuations occur at q = n corresponding to the 
staggered current pattern of the staggered flux state. The 
same behavior is also observed for higher doping as can 
be seen in Fig. || showing similar data but for a doping 
of four holes. 

The behavior of the integrated spectral weight, namely 
the current-current susceptibility, as a function of param- 
eters and doping is also of particular interest. Data re- 
ported in Fig. H indeed reveal that a qualitative change 
also occurs at the transition between bound hole pairs 
and free holes. When holes are tighly bound into pairs 
(typically for t' > 0) the current susceptibility is rather 
featureless and scales with the number of holes. When 
t' — » t' c from above, broad peaks appear in the vicinity 
of momentum ir flffl especially in the system containing 
only two holes. This increase of total spectral weight 
is mainly due to the increase of the low energy incom- 
mensurate peaks in the dynamical correlation functions 
shown in Figs. §(b), ||(b). When hole pairs break apart 
i.e. for t' < t' c a large peak appears at the commensurate 
wavevector q = n. 

To complete our study, we have also computed 
the equal-time current-current correlation function 
J2 q e%qr '{jqj-q) vs separation r between two rungs. Re- 
cent Density Matrix Renormalisation Group (DMRG) 
studies on a 2 leg t-J ladder have shown that such corre- 
lations are exponentially suppressed with distance in the 
LE liquid phase . Our data for various t' and 2 differ- 
ent hole concentrations are plotted in Fig. I| Above t' c , in 



FIG. 3. Dynamical correlation function Sjj(q,uj) calcu- 
lated on a 4 hole-doped periodic ladder for J = 0.5. Same 
presentation as Fig. ^ 

the LE liquid phase, the alternating current correlations 
are indeed rapidly decreasing with distance in agreement 
with the DMRG results. When holes are strongly bound 
into pairs (as it is the case e.g. for t' = 0.4) it is inter- 
esting to notice that, apart from an overal scaling factor 
proportional to the hole concentration, the current cor- 
relation is almost completely independent of the number 
of pairs in the system. This is a further confirmation 
that current fluctuations occur within one pair. On the 
contrary, when pairs break apart the current correlations 
within the GS increase strongly |L9| with decreasing hole 
density. In this case, establishment of staggered current 
correlations depends strongly on the parity of the ladder 
length because of possible frustration as seen in Fig. ||(a) . 

We finish this paper by briefly discussing the role of 
the SU(2) symmetry of the model. As mentioned above, 
it leads to the equivalence between d-wave pairing and 
staggered flux phases at half-filling [Q. Since this de- 
generacy is lifted upon doping, the hole kinetic energy 
is expected to stabilize one phase over the other. The 
sudden change of properties in a t-t'-J ladder at small 
doping when varying t' is consistent with this scenario 
if one assumes that the competition between the pairing 
and the staggered flux phases is indeed governed by the 
small kinetic energy term (proportional to doping) and 
hence by t' . 

In summary, our data strongly suggest that, when hole 
pairs are stable, the current fluctuations is intrinsically 
related to the internal spin/charge dynamics within one 
pair at least in the regime where pairs do not overlap. 
The magnitude of the current fluctuations increases as 
the spatial extension of the pairs increases until they start 
to overlap. These findings are consistent with the picture 
of bound hole pairs of holes carrying opposite current 
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FIG. 4. Current-current susceptibility vs momentum q for 
different NNN hopping t' (as indicated on the plots) calcu- 
lated on 2 and 4- hole doped 2x 12 ladders. The susceptibilities 
have been normalised w.r.t. the number of holes to emphasize 
proportionality if any. Arbitrary units are used. Note that at 
t' = 0.4 data for 2 and 4 holes are almost indistinguishable. 



FIG. 5. Current (equal-time) correlation function vs dis- 
tance for two (a) and four (b) hole-doped ladders calculated 
on a 2x12 cluster. Different filled symbols correspond to dif- 
ferent t' parameters according to the legend. In (a) data for 
a 2x11 ladder are also shown as open symbols (unless undis- 
tinguishable from the previous sets). 



vorticity §. When hole pairs break apart into individ- 
ual holes, commensurate staggered fluctuations set in. 
The increase of the magnitude of these fluctuations as 
one moves towards vanishing doping is interpreted as an 
instability to a (static) staggered flux state. 
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